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(g) Method and apparatus for measuring concentration of absorptive constituent in scattering medium. 



@ An optical pulse emitted from an optical pulse 
source Is incident on a living tissue and is input 
to an optical sampling unit The optical sampl- 
ing unit detects intensities f of optical pulses 
detected at the respective timings of trigger 
signals Input from a delay unit A first ao- 
cumulator calculates products t f of the inten- 
sities f and the delay times t from the delay unit 
A second accumulator accumulates the light 
intensities f. An average optical pathlength cal- 
culating unit calculates an average optical 
pathlength using the accumulation results as 
the products t f and the intensities f. An SOj 
value calculating unit calculates a ratio of VHb02 
to from the average optical pathlength to 
calculate an SO2 value. 
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BACKGROUND OF THE INVENTION 
Field of The Invention 

5 The present invention relates to a method and apparatus for measuring the concentration of an absorptive 

constituent (e.g., hemoglobin) present in a scattering medium (e.g., a living tissue). 

Related Background Art 

10 When light having a wavelength near a near-infrared range is incident on a living tissue, absorption of the 

incident light In the living tissue is greatly influenced by hemoglobin in blood. The influence of oxyhemoglobin 
(HbOa) is different from that of deoxygenated hemoglobin (Hb). The absorption spectra of the living tissue 
change in accordance with the states of hemoglobin, as shown in a graph of FIG. 1. The wavelength [nm] of 
incident light is plotted along the abscissa in FIG. 1, and the absorbance [mM-icm-^] is plotted along the or- 
is dinate therein. As shown in this graph, the absorbance caused by deoxygenated hemoglobin (Hb) is higher 
than that by oxyhemoglobin (HbOj) when the wavelength is shorter than 800 nm. When the wavelength is lon- 
ger than 800 nm, the absorbance caused by oxyhemoglobin (HbOj) is higher than that by deoxygenated he- 
moglobin (Hb). That is, the optical absorption in the living tissue is understood to change in accordance with 
the metabolic conditions of oxygen in the living tissue. 

20 Conventionally, a monitor for oxygen metabolism in a living tissue has been developed using this principle. 

Such a monitor has been put to practical use. In this oxygen metabolism monitor, however, a relative change 
from the start of measurement can only be monitored. For this reason, the degree (SO2 value) of saturation 
of oxygen in blood in a living tissue, which is a ratio of an oxyhemoglobin (Hb02) concentration to a total her 
moglobin concentration, cannot be measured. Therefore, the SO2 value as the index of oxygen metabolism in 

25 the living tissue cannot be known. Therefore, an attempt has been made using the following method and ap- 
paratus to measure this SO2 value. 

First, a measuring apparatus (tradename: Cerebral Oximeter (invos 3100)) available from Somanetics is 
used as an apparatus for measuring an SO2 value in blood in a living tissue. This apparatus nrieasures the 
concentration of an absorptive constituent (hemoglobin) in the scattering medium (living tissue) using the prin- 

30 ciple shown In Fig. 2. More specifically, continuous (CW) light having different wavelengths in the near-infrared 
range is incident on a living tissue 1, and attenuation amounts of the incident light are detected at detection 
positions spaced apart from the light incident position by distances r1 and r2. The Hb and HbOa concentrations 
are obtained In accordance with the conrelations between the attenuation amounts of incident light components 
having different wavelengths and the distances r1 and r2 between the light incident position and the detection 

35 positions, because the absorption profiles of incident light components in the living tissue is different as de- 
scribed above, thereby measuring the SO2 value in blood. 

Second, a spectroscopic technique for an absorptive constituent in a scattering medium using a time- 
resolved spectroscopy Is available, as disclosed in U.S. Patent No. 5,11 9,81 5 by B. Chance et al. This reference 
describes that the spectroscopic technique is applied to an SO2 value measurement in blood of a living tissue. 

40 More specifically, an optical pulse is incident on the living tissue, and an optical pulse profile spreading as a 
function of time due to light scattering is time-resolved measured to obtain a profile representing a change in 
light intensity as a function of a change in time. Ught absorption in the Irving tissue is measured such that the 
light intensity of the resultant profile is logarithmically calculated, and a gradient of the change in light intensity 
as a function of time is obtained. When time-resolved measurement is performed upon incidence of two optical 

45 pulses having different wavelengths on the living tissue, the light absorption for each optical pulse of each wa- 
velength is measured, thereby the SO2 value in the blood can be calculated, because the light absorption pro- 
files of Hb and Hb02 are different from each other. 

Third, a spectroscopic technique for an absorptive constituent In a scattering medium, using phase mod- 
ulation spectroscopy, is available, as disclosed in U.S. Patent No. 4,972,331 by B. Chance et al. This reference 

50 describes that this spectroscopic technique can be applied to an SO2 value nrieasurement In blood. More spe- 
cifically, modulated light is incident on the living tissue, and light absorption of the incident light Is detected on 
the basis of a change In phase caused by propagatton of the modulated light through the living tissue. When 
modulated light components having different wavelengths are Incident on the living tissue, the concentration 
ratios of Hb02 and Hb are detected to calculate the SO2 value in the blood because the change in phase varies 

55 depending on the type of absorptive constituent in the living tissue and the wavelength of the incident modu- 
lated light 
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SUMMARY OF THE INVENTION 

The following problems are posed by the conventional measurements for concentrations of absorptive con- 
stituents In the scattering media. 

The measuring apparatus available from Somanetlcs, as the first prior art. must have a plurality of pho- 
todetecting devices to result in complicated device setup because the attenuation amounts of the incident light 
must be detected at a plurality of detection positions. In addition, conditions under which the respective pho- 
lodetecting devices are brought into contact with a skin must be equal to each other. When a patient moves 
in contact with the photodetecting devices even under the equal conditions, photodetecting conditions of the 
respective devices may become nonuniform at the respective positions. That is, indefinite factors are included 
in optical measurements in this measuring device. It is, therefore, difficult to accurately measure ah SO2 value 
in blood. 

In use of the time-resolved measurement as the second prior art, the total light intensity vs. time profile 
must be measured because the optical absorption of the incident light is detected in accordance with the gra- 
dient of the light intensity vs. time profile. For this reason, the calculation process for obtaining the SO2 value 
from the measurement data becomes complicated, and a circuit arrangement required for this calculation proc- 
ess becomes bulky. In this measurement method, the number of wavelengths of the optical pulses Incident on 
the living tissue is two. arid the incident pulse light having each wavelength is time-resolved measured to cal- 
culate the SO2 value in the blood. Such a two-wavelength measurement can be performed only when back- 
ground absorption except for hemoglobin in the blood can be neglected. H2O, proteins, and the like present 
in the living tissue except for the hemoglobin cannot be generally neglected with respect to incident light In the 
near-Infrared range. Therefore, the concentration measurement of the absorptive constituents in the scattering 
medium, proposed by B. Chance et al.. can be applied to only a limited number of scattering media so as tO: 
measure SO2 values. 

In the technique for measuring an SO2 value using the frequency-resolved measurement method as the 
third prior art It is expected that an optimal modulation frequency of the incident light may vary depending on 
scattering condition, boundary conditions, and the like In the living tissue. For this reason, the modulation fre- 
quency must be set again for each living tissue on which an optical pulse is incident. The optical measurement 
process cannot be simply performed. In additbn, when the modulation frequency cannot be optimally set for 
the living tissue, an accurate SO2 value cannot be measured. Even in this prior art using the modulated light, 
a two-wavelength infrared measurement is required. The background absorption cannot be neglected, as de- 
scribed above, and this technique can be applied to only a limited number of scattering media so as to measure 
the concentrations of absorptive constituents. 

The present invention has been made to solve the conventional problems described above and is charac- 
terized by comprising causing light having a predetermined wavelength to be incident from a predetermined 
position of a scattering medium containing n (n ^ 2) types of absorptive constituents, detecting the light, scat- 
tering in the scattering medium, at a predetermined position spaced apart from the light incident position by a 
predetermined distance, calculating an average optical pathlength of the light having the predetermined wa- 
velength in the scattering medium on the basis of the detected light, obtaining respective average optical path- 
lengths of light components having n+1 types of wavelengths while a wavelength of light incident on the scat- 
tering medium is changed, obtaining respective average optical pathlength differences between respective 
pairs of light components having different wavelengths from the average optical path lengths, and obtaining 
a concentration ratio of the respective absorptive constituents in the scattering medium In accordance with a 
relationship representing that the average optical pathlength difference is in inverse proportion to a difference 
In absorbance of the absorptive component between the light components having two wavelength. 

The present invention is also characterized by comprising a light source for generating pulse light having 
a predetermined wavelength every predetermined time interval, light guide means for causing each pulse light 
from the light source to be incident on a scattering medium containing n (n s 2) types of absorptive constituents 
from a predetermined position and extracting each pulse light scattering within the scattering medium from a 
predetermined position spaced apart from the light Incident position by a predetermined distance, photode- 
tecting means for tlme-resdved measuring each pulse light extracted by the light guide means, wavelength 
control means for controlling a wavelength of the pulse light emitted from the light source, optical pathlength 
calculating means for calculating an average optical pathlength of light components having n+1 wavelengths 
in the scattering medium in accordance with time-resolved measurement results of pulse light components 
having n+1 wavelengths obtained by causing the wavelength control means to control the wavelength of the 
pulse light emitted from the light source, and concentration calculating means for obtaining an average optical 
pathlength difference between light components having two wavelengths from the average pathlengths cal- 
culated by the optical pathlength calculating means and obtaining a ratio of the absorptive constituent con- 
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centrations in the scattering medium in accordance with a relationship representing that the average optical 
pathlength difference is in Inverse proportion to a difference In absorbance of the absorptive component be- 
tween the light components having two wavelengths. 

According to the present invention, incident light scattering in a scattering medium can be detected at one 
5 location of the scattering medium. Therefore, detection of each incident light passing through the scattering 
medium can always be performed under the same conditions, and indefinite factors in the optica) measurement 
process can be excluded. For this reason, the concentration measurement in the scattering medium can be 
accurately performed. 

Even If background factors except for absorptive constituents serving as measurement targets are present 
10 in the scattering medium, the measurement is free from the influences of light absorption caused by these 
background factors. Therefore, the concentration measurement In a wide range of objects can be performed 
regardless of the types of scattering media. 

Light incident on the scattering medium can be set independently of light scattering and the like in the scat- 
tering medium. For this reason, the conditions of the apparatus need not to be changed in accordance with 
15 the types of scattering media, and hence the conditions of the incident tight need not to be set again. The con- 
centration measurement of the absorptive constituent in the scattering medium can be easily and accurately 
performed. 

BRIEF DESCRIPTION OF THE DRAWINGS 

20 

Fig. 1 is a graph showing the light absorption profiles of deoxygenated hemoglobin Hb and oxyhemoglobin 
HbOz. 

Fig. 2 Is a sectional view showing a living tissue for explaining the conventional principle of measuring the , 
concentration of an absorptive constituent in a scattering medium. 
25 Fig. 3 shows that pulse light scatters In a living tissue spreading as a function of time and is detected. 

Fig. 4 shows a graph for explaining the principle of measurement of the present invention, representing a 
profile X of pulse light incident on a scattering medium and profiles A and B of pulse light components detected 
upon scattering and absorption of this pulse light in the scattering medium. 

Fig. 5 is a graph showing sampling timings t„ of a detected light profile Y when an average optical path- 
30 length is expressed in the form of a sequence. 

Fig. 6 is a graph showing a detected light profile Y obtained when a large number of photons are detected 
in single-photon measurements. 

Fig. 7 is a block diagram showing the schematic arrangement of a light Intensity sampling apparatus ac- 
cording to the f rst embodiment of the present invention. 
35 Fig. 8 is a timing chart showing signals In the respective components of the light Intensity sampling appa- 
ratus according to the first embodiment. 

Fig. 9 is a block diagram showing the schematic arrangement of a light intensity sampling apparatus ac- 
cording to the second embodiment of the present invention. 

Fig. 10 is a timing chart showing signals in the respective components of the light intensity sampling ap- 
40 paratus according to the second embodiment 

Fig. 11 is a block diagram showing the schematic arrangement of a light intensity sampling apparatus ac^ 
cording to the third embodiment of the present invention. 

Fig. 12 Is a view showing a rat serving as a sample used in the third embodiment. 
Fig. 13 Is a graph showing the profiles of incident pulse light and detected pulse light in the third embodi- 
45 ment. 

Fig. 14 is a graph showing changes in average optical pathlength with changes In wavelength of incident 
pulse light when an FIO2 concentration In an inhalatbn gas supplied to the rat Is used as a parameter. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

50 

A method and apparatus for measuring the concentration of an absorptive constituent In a scattering me- 
dium according to the present Invention will be described when they are applied to a measurement of Hb and 
Hb02 concentrations in a blood of a living tissue. In principle, this concentration measurement comprises [1] 
a process for correlating an average optical pathlength in the scattering medium with light absorption informa- 
55 tion. [2] a process for measuring an average optical pathlength, [3] a process for acquiring light absorptfon In- 
formation from an average optical pathlength difference between optical pulses having two different wave- 
lengths, and [4] a process for deriving a ratio of absorptive constituent concentrations in the scattering medium 
from this light absorption information. These processes will be described In detail below. 
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[1] Process for Correlating Average Optical Pathlength with Light Absorption Information 

As shown in Fig. 3, when an optical pulse is incident on a living tissue 11, the light pulse scatters as a 
function of time and is detected. The intensity of the optical pulse detected as described above has a profile 

5 shown in a graph in Fig. 4. Time [ns] is plotted along the abscissa in Fig. 4, and the detection light intensity 
[a.u.] is plotted along the ordinate therein. As shown in this graph, the average optical pathlength of the optical 
pulse between the optical pulse incident position and a detection position changes in accordance with an in- 
fluence of light absorption in the living tissue. More specifically, when pulse light having a profile X indicated 
by a doted line In Fig. 4 Is incident on the living tissue, a detection prof ile is detected as a profile A when light 

io absorption In the living tissue is high; when light absorption is low, the detection profile becomes a profile B. 
That is, a mean flight time <t> during which the optical pulse passes through the living tissue changes in ac- 
cordance with a light absorption state in the living tissue. This mean flight time <t> corresponds to a time in- 
terval from an incident timing at zero time when the pulse light is incident to a barycentral position a or b of 
the profile A or B. As shown in this graph, when the light absorption is high, the mean flight time <t> becomes 

15 short; when the light absorption is low, the mean flight time <t> becomes long. The present inventors disclose 
in other patent applications (Japanese Patent Application No. 2-322105 and Japanese Patent Laid-Open No. 
4-1 91 642) that a detection intensity f(t,8\/C) actually measured for an incident optical pulse having such spread- 
ing as a function of time is represented as a product of a light intensity fo(t) in the absence of absorptive con- 
stituents in the scattering medium and a light Intensity iO-^vct jn the presence of absorptive constituents as 

20 follows: 

f(t,eVC) = fo(t)10- «va (i) 

where 

t is the time from the optical pulse Incident timing 

E is the absorption coefficient of an absorptive constituent in the scattering medium, which varies de- 
25 pending on wavelengths of incident light components 

V is the molar concentration of the absorptive constituent 
C is the speed of light in the scattering medium 

Equation (1 ) is deemed to represent a photon distribution as a function of time until the photon reaches 
from the light incident position to the detection position. That is, this profile is a so-called scattering time profile. 
30 The parameters £, V, and C inaease with respect to the scattering time profile fo(t) in the absence of light ab- 
sorption. To siniplify the calculation, the exponential base of 10 is replaced with e, and the exponential part - 
sVCt is replaced with -at to rewrite equation (1) into equation (2) as follows: 

f(t,a) = fo(t)e - «t (2) 

In this case, the absorption coefficient for the exponential base of e is defined as e*. and eVC = 0.434eVC 
35 = a is established. 

An average optical pathlength C<t(a)> required to cause the photon to reach the detection position from 
the light Incident position is defined as equation (3) below: 

^ cfcf^(c)e'^'dC 

C<C(a)>=— ^ ...(3) 

45 * 3 



50 



The numerator and denominator of equation (3) are Taylor-expanded to approximate terms up to terms of 
the first degree, and the resultant equation is divided by an integral term represented by following expression: 



fAOe'^^'dt ...(4) 



55 Equation (5) can be obtained by the following equation transfomn. In this case, a Laplasian transform re- 

lation is used to derive equation (5), provided that a « do + Aa is established: 
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10 



IS 



<C{a)> = 



_ 0 



0 

0 

0 . _o 

0 0 

<c(<3o)> - (a-aa)<c-(a-,)> 



(5) 



20 ••• <t(a)> 



1 - (a-a,) <t(ao) > 

<f:(aa) > - <f:'(ao)>Aa 
1 - <t(aQ)>Aa 



[2] Process for Measuring Average Optical Pathlength 

25 The average optical pathlength is defined by equation (3), but can be rewritten as equation (6) in the form 
of a sequence. 

n 

wherein times ti. t2. ta, ...) can be plotted at a predetermined interval, as indicated in the graph of 
Fig. 5. Time t is plotted along the abscissa of this graph, and a light intensity f(t) is plotted along the ordinate 
35 thereof. This graph represents that an optical pulse having a profile X indicated by a dotted line is incident on 
a living tissue to obtain a detection profile Y. 

The average optical pathlength C<t> is obtained from equation (6) as follows. That is, timings at which 
optical pulses are incident on the living tissue are given with delay times tn with reference to the incidence of 
the light having the profile X, and light intensities f(tn.eCV) are measured at the respective timings. The products 
40 of the light intensities f(tn,eCV) and delay times tn, which are represented as the numerator of equation (6), are 
accumulated. The light intensities f(tn,eCV) represented by the denominator of equation (6) are accumulated. 
The sum of the numerator is divided by the sum of the denominator to obtain a quotient, thereby calculating 
the average optical path length C<t>. 

The average optical pathlength can also be obtained in a single-photon measurement by equation (7) as 
45 follows: 



C<C(a) > = 



so 



M 



. .(7) 



where N is the number of measured photons. 
Adetection profile shown in Fig. 6 is obtained according to this single-photon measurement As in the graph 
of Fig. 5, time t is plotted along the abscissa, and the light intensity f(t) is plotted along the ordinate in Fig. 6. 
55 When an optical pulse having a profile X is Incident on a living tissue, and a plurality of photons having different 
delay times from the optical pulse incident timing are respectively measured, a detection profile Y is obtained. 
The delay times from the incident timings of individual photons to their detection timings are accumulated, and 
the sum is divided by the total number N of photons, thereby calculating the average optical pathlength. 
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[3] PROCESS FOR ACQUIRING LIGHT ABSORPTION INFORMATION FROM AVERAGE OPTICAL PATH- 
LENGTH DIFFERENCE 

Light absorption information In the scattering medium is obtained from the resultant average optical path- 
5 length as follows. That Is, equation (5) is transformed into equation (8) below: 

<t(a)> - <t(a)><t(ao)>Aa = <t(ao)> - <t2(ao)>Aa 
/. <t(a)> . <t(ao)> = (<t(a)><t(ao)> - <t2{ao)>)Aa (8) 
In this case, when an approximation is performed such that <t(a)><t(ao)> = <t(ao)>2, equation (8) can be 
rewritten as equation (9) as follows: 

<t(a)> - <t(ao)> = (<t(ao)>2 - <t2(ao)>)Aa 
/. <t(a)> - <t(ao)> = P(ao)Aa (9) 
for P(ao) = <t(ao)>2 - <t2(ao)>. wherein this P(ao) is negative. 
When optical pulses having two different wavelengths X1 and X2 and equal scattering coefficients in the 
scattering medium, i.e., equal values fo(t). are incident on the living tissue, equation (10) can be derived from 
IS equation (9). In this case, the variable a at the wavelength XI , the constant ao at the wavelength , the variable 
a at the wavelength X2, and the constant ao at the wavelength X2 are defined as Bxa, Qoxu a^a. and aox2. 

. <t(axi)> - <t{aoxi)> = P(aoxi)Aaxi 
<t{ax2)> - <t(ao«)> = P(aox2)Aa« (10) 
In equations (1 0). aoxi and ao^ depend on scattering. If aox^ = aou = a^ then Aa^i = a^i- ao and Aax2 = b^z 
20 - ao, thereby obtaining equation (11 ) as follows: 

<t(a;^> - <t(a,,)> = P(ao)- (a« - a„) (11) 
Since eVC = a is established, a^i = sxiCV and ax2 = exaCV where e^i is the absorption coefficient at the 
wavelength X1, and is the absorption coefficient at the wavelength X2. When the exponential base is re; 
turned from e to 10. equation (11) can be rewritten as equation (12). provided that P(ao) with the exponential 
25 base of 10 is represented as P'(ao): 

{<t(^CV)> - <t(exiCV)>} = P'(ao)(ex2 - eii)VC (12) 
This equation (12) represents that the difference between the average optical pathlengths or the mean 
flight times of the optical pulses having two different wavelengths is in Inverse proportion to the absorbance 
in the scattering medium. Therefore, the light absorption Information in the scattering medium can be obtained 
30 by calculating the average optical pathlength difference or the mean flight time difference. 

[4] PROCESS FOR DERIVING RATIO OF ABSORPTIVE CONSTITUENT CONCENTRATIONS IN SCAT- 
TERING MEDIUM 

35 A concentration ratio of the absorptive constituents in the scattering medium, i.e., the SO2 value in the 
blood of the living tissue, is obtained from the average optical pathlength difference calculated as descrit)ed 
above. For this purpose, a concentration ratio of oxyhemoglobin (HbOz) to deoxygenated hemoglobin (Hb) 
must be obtained. Two average optical pathlength differences represented by equation (12) must be calculat- 
ed, and measuremeints of the average optical pathlengths at least three wavelengths (X1, X2, and X3) must 

40 be perfomied. The absorption coefficients of HbOa and Hb at the respective wavelengths are summarized in 
Table 1 below. 



[Table 1] 



45 



Wavelength 


Hb02 


Hb 


XI 




€Hb.X1 


12 






X3 




eHb,X3 



Assume that the HbOj and Hb concentrations In the living tissue, and the mean flight times during which 
the optical pulses at the wavelengths XI , X2, and X3 pass through the living tissue, defined as VHb02. VHb. <t>:xi , 
<t>x2, and <t>x3, respectively. Light absorption is also present due to background factors except for hemoglobin 
55 in the living tissue. Assume light absorption values caused by these background factors at the wavelengths 
XI, X2, and X3 are defined as a^i, a^a, and a^a, respectively. At this time, light absorption a^i, a,^, and a^a upon 
incidence of the pulse light components having the wavelengths XI , X2, and X3 are defined by equations (1 3) 
below: 
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3X1 = (eHb02A1 * VHb02 ^HbAI ' ^Hb + 
3X2 = (£Hb02,X2 ' VHb02 ^ eHb,X2 " ^Hb + ^iJC 
3X3 = (£Hb02.X3 • VHb02 + eHb,X3 ' - aj^)C (13) 

Equation (12) can be rewritten as equation (14) below. In this case, assume the background absorption 
5 values in the living tissue at the wavelengths XI , A2, and X3 are equal to each other (a^i = = axa). 

<t>x2 - <t>xi = P'(ao){(eHb02.K2 - eHb02.Xl)VHb02 + (£Hb,X2 - eHb,xi)VHb}C 
<t>W - <t>Xl = P'(ao){(€Hb02,X3 - eHb02.Xl)VHb02 + (£Hb^ - eHb.Xl)VHb}C (14) 

In this case, if <t>x2 - <t>xi = A<t>x2,xi and <t>x3 - <t>xi = A<t>x3,xi, equations (1 5) can be derived as follows: 

(eHb02,X2 - eHb02,x.i)VHb02 + (eHb.X2 - eHb,xi)VHb = A<t>x2.xi/P*(3o)C 

(eHb02.X3 - £Hb02.Xt)VHb02 + (eHb.X3 " ^Ht>^^)'^Ht> = A<t>x3.Xl/P'(ao)C (15) 

Equations (15) are solved as simultaneous equations using VHboz and as unknown values to obtain 
the following equations (1 6):. 

PXao)VHbo2 = {(eHb.x3 - eHb.xi)A<t>x2.xi " (eHb,x2 " ^b;K^)^Ow.^)fO 

P'{ao)VHb = {-(eHb02.X3 - eHb02,Xl)A<t>X2A1 + (eHb02Jl2 - eHb02.3Cl)A<t>x3.Xl)/D (16) 

15 In this case, D is represented as follows: 

D = {(eHb02.X2 - eHb02.Xl)(€Hh;c3 - CHbAl) - (fiHb^ " SHbAl)(fiHb02.X3 " eHb02Al)}C (^7) 

The ratio of the oxyhemoglobin concentration VHb02 to the deoxyhemoglobin Vh, can be obtained from 
equations (1 6) as follows: 

VHb02:VHb= {{SHhXA " €Hb.3ll)A<t>X2,X1 - (CHb,X2 ' eHb.Xl)A<t>x3.xi} : { -(eHb02Jl3 - €Hb02Al)AOx2^1 + (eHb02^ " 

2^ eHb02,xi)A<t>x3.xi} (18) 

The degree (SO2 value) of saturation of oxygen is calculated as follows: 

SO2 = VHb02/(VHb02 + VHb) = 1/(1 + VHb/VHb02) (19) 

A ratio of VHb02to VHb is obtained from equation (1 8). and a substitution of this concentration ratio in equa-. - 
tion (19) can yield the degree of saturation of oxygen. 
25 The preferred embodiments of the present invention using the above principle of measurement will be de- 
scribed below. 

Fig. 7 is a block diagram showing the schematic arrangement of a light intensity sampling apparatus ac- 
cording to the first embodiment of the present invention. This light intensity sampling apparatus measures he- 
moglobin concentrations in the living tissue. This apparatus will be described with reference to the timing chart 
30 in Fig. 8. 

An optical pulse source 21 emits an optical pulse having a predetermined wavelength XI at a predeter- 
mined time interval (see a part (a) of Fig. 8). The emitted optical pulse is incident on a living tissue 23 through 
an optical fiber 22 or directly. The optical pulse spreading as a function of time due to light scattering in the 
living tissue 23 is Input to an optical sampling unit 25 directly or through an optfcal fiber 24. The optical pulse 

35 source 21 outputs a trigger signal to a delay unit 26 at an optical pulse emission timing. The delay unit 26 delays 
the input trigger signal by a predetermined period of time and outputs the delayed trigger signal to the optical 
sampling unit 25. This delay time is appropriately changed in each optical pulse. That is, the delay time changes 
in an order of ti, t2, tg, t4, ... with reference to the optical pulse emission timings (see a part (b) of Fig. 8), 
The optical sampling unit 25 receives each delayed trigger signal from the delay unit 26. At the respective 

40 sampling timings when the connesponding delayed trigger signals are input, the optical sampling unit 25 detects 
Intensities fi. fa, fa. U*— of the optical pulses emitted from the optical pulse source 21 and passing through the 
living tissue 23 (see a part (c) of Fig.8). The delay unit 26 outputs the delay time of each delayed trigger signal 
to a first accumulator 27. The first accumulator 27 calculates a product t f of the detected light intensity f sup- 
plied from the optical sampling unit 25 and the delay time t supplied from the delay unit 26 at each sampling 

45 timing. Such products tf are accumulated (trfi + tj fa + tg-fg + t4 f4. ...). and the sum is output to an average 
optical pathlength calculating unit 29 (see a part (d) of Fig. 8). A second accumulator 28 receives the intensities 
f of the optical pulses, detected through the living tissue 23, from the optical sampling unit 25. The received 
intensities f are accumulated at the respective sampling timings (fi + fa + fa + U, ...), and the sum Is output to 
the average optical pathlength calculating unit 29. (see a part (e) of Fig. 8). 

50 The average optical pathlength calculating unit 29 calculates an average optical pathlength based on equa- 
tion (6) using the sum of the products tf Input from the first accumulator 27 and the sum of the light Intensities 
f Input from the second accumulator 28. This average optical pathlength is obtained from an optical pulse 
spreaded and detected as a function of time upon incidence of the optical pulse having the predetermined wa- 
velength XI on the living tissue 23. The average optical pathlength calculating unit 29 outputs a control signal 

55 to the optical pulse source 21 to change the wavelength of an optical pulse emitted from the optical pulse 
source 21 . By this wavelength control, the above calculation processes are also performed for the wavelengths 
X2 and X3 of the optical pulses, and the average optical pathlength calculating unit 29 also calculates average 
optical pathlengths for these wavelengths X2 and X3. 
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An SO2 value calculating unit 30 receives the average optical pathlength calculation results for the wave- 
lengths X.1 , X2, and X3 and calculates a ratio of the oxyhemoglobin concentration VHb02 to the deoxyhemoglobin 
concentration Vnb in the blood of the living tissue on the basts of equation (18). The resultant concentration 
ratio is substituted into equation (19) to calculate an SO2 value in the blood. This SO2 value is displayed on a 
5 display unit 31. 

According to this embodiment, the pulse light emitted from the optical pulse source 21 and scattering in 
the living tissue 23 is detected at one location of the living tissue 23. Unlike the conventional measuring ap- 
paratus available from Somanetics. incident pulse light need not to be detected at the plurality of positions of 
the scattering medium. For this reason, once the pulse light incident position and the incident pulse light de- 
10 tection position are set, the respective pulse light components can always be detected under the same con- 
ditions. The measurement is free from influences of a contact state of a photodetecting unit with a skin and a 
motion of the living tissue. Therefore, the concentration measurements can always be accurately performed. 

Even if background factors such as H2O and proteins except for hemoglobin are present, the measurement 
is free from the influences of light absorption caused by these factors. A ratio of Hb02 to Hb in the living tissue 
IS can be obtained without receiving these Influences. For this reason, according to this embodiment, the con- 
centrations of absorptive constituents in a wide range can be measured without Influencing the types of scat- 
tering means. 

Pulse light incident on the living tissue 23 can be set independently of light scattering in the scattering 
medium. Unlike the prior art disclosed in U.S. Patent No. 4,972,331, the apparatus need not to be reset de- 

20 pending on the types of scattering media as measurement targets or the incident pulse light conditions need 
not to be changed. For this reason, the concentration measurement of the absorptive constituent in the living 
tissue can be facilitated, and a problem posed by inappropriate setup of the apparatus to lead to an inaccurate 
concentration measure can be prevented. 

A concentration measurement of absorptive constituents in a scattering medium according to the second 

25 embodiment of the present invention will be described below. Fig. 9 is a block diagram showing the schematic 
arrangement of a light intensity sampling unit according to the second embodiment, and this unit will be de- 
scribed with reference to Fig. 10. 

An optical pulse source 41 emits an optical pulse having a predetermined wavelength X^ at a predeter- 
mined time interval (see a part (a) of Fig. 10). The emitted light pulse is incident on a living tissue 43 through 

30 an optical fiber 42 or directly. The light pulse propagating through the living tissue 43 is Input to an photomul- 
tiplier tube (PMT) 45 through an optical fiber 44 or directly. When the optical pulse emitted from the optical 
pulse source 41 is input to this PMT 45. one or less photon must be detected per optical pulse single photon 
event. This one photon is converted into an electrical pulse signal by the PMT 45. 

This electrical pulse signal is input to a constant fraction discriminator (CFD) 46 and the timing signal is 

35 picked off. The photons detected by the PMT 45 have different delay times from the optical pulse emission 
timings for the respecth/e optical pulses. The CFD 46 outputs pulse signal (see a part (b) of Fig. 10) delayed 
from the optical pulse emission timing by times T,, T2, T3, T4, .... Note that a wavefomn indicated by a dotted 
line in the part (b) of Fig. 10 represents the profile of detection optical pulses obtained upon detection of a large 
number of photons. This profile represents spreading of the input optical pulse as a function of time. 

40 The output pulse form the CFD 46 is supplied to a time to amplitude converter (TAC) 47. A timing signal 
is supplied to the TAC 47 when the optical pulse is emitted from the optical pulse source 41 . The TAC 47 con- 
verts time differences between the timings of optical pulse incidence on the living tissue 43 and the detection 
timings of the optical pulses in the PMT 45 on the basis of the pulse signals from the CFD 46 and the timing 
signals from the optical pulse source 41. That is. the TAC 47 outputs a pulse signal having an amplitude cor- 

45 responding to each time difference (see a part (c) of Fig. 10). 

A sample & hold (S & H) circuit 48 receives an output pulse from the TAC 47 to detect and hold the peak 
value of this signal. An analog-to-digltal converter (ADC) 49 converts this peak value into a digital signal. This 
value con*esponds to each of the amplitudes Ti to T4 (see a part (C) of Fig. 10) of the electrical pulse signals. 
An accumulator 50 accumulates these digital values (T, + T2 + T3 + T4, ...) to obtain the sum of the amplitude 

50 values of the respective electrical pulse signals (see a part (d) of Fig. 10). The accumulator 50 counts the nunv 
ber of electrical pulse signals output from the TAC 47 In addition to the above accumulatton (see a part (e) of 
Fig. 10). These accumulation values are output to an average optical pathlength calculating unit 51. 

The average optical pathlength calculating unit 51 divides the sum of the amplitudes of the respective elec- 
trical pulse signals, i.e., the sum total of the photon detection delay times, by the total photon count to obtain 

55 a quotient, and outputs this quotient as an average optical pathlength to an SO2 calculating unit 52. The aver- 
age optical pathlength calculating unit 51 outputs a control signal to the optical pulse source 41 to change the 
wavelength of each optical pulse emitted from the optical pulse source 41. The above calculation processes 
are also performed for the wavelengths X2 and 7J3 of the optical pulses, and the average optical pathlength 
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calculating unit 51 also calculates the average optical pathlengths for the wavelengths X2 and X3. The SO2 
calculating unit 52 receives the average optical pathlength calculation results of the wavelengths XA, X2, and 
X3 and calculates a ratio of an oxyhemoglobin concentration VHb02 to a deoxyhemoglobin concentration 
in the blood of the living tissue, on the basis of equation (18). The resultant concentration ratio is substituted 
into equation (1 9) to obtain an SO2 value. The SO2 value is output to a display unit 53. The display unit 53 dis- 
plays the input SO2 value. 

The same effect as in the first embodiment can be obtained in the second embodiment The degree of 
saturation of oxygen in the blood of the living tissue can be easily and accurately measured with a simple ap- 
paratus arrangement. 

The third embodiment will be described below in which a concentration measurement of absorptive con- 
stituents in a scattering medium Is applied to a rat. Fig. 11 is a schematic block diagram of a concentration 
measuring apparatus according to the third embodiment 

A mode-locked Ti:Sapphire laser 61 is a pIco-second optical pulse source and is excited by an output from 
an Ai^ laser 62. In this embodiment, laser pulse light (pulse width: 2 pS or less) having a repetition frequency 
of 76 MHz within the wavelength range of 730 to 860 nm is output from the mode-locked Ti:Sapphire laser 61 . 
This pulse light is incident on a sample 64 through an optical fiber 63. The pulse light diffusively reflected by 
or passing through the sample 64 is input to a streak camera 66 through an optical fiber 65. At the same time, 
the pulse light output from the Ti:Sapphire laser 61 is input to the streak camera 66 through the optical fiber 
67 without going through the sample 64. In this embodiment, the pulse light incident on the sample 64 is time- 
resolved measured by the streak camera 66. 

The sample 64 is a Wislar rat shown in Fig. 12. The pulse light incident optical fiber 63 is inserted into the 
oral cavity of the rat 64. The pulse light detection optical fiber 65 Is located at the parietal portion of the rat 64. 
This rat 64 is anesthetized with nembutal, and a tracheal cannula is inserted into a trachea. An artificial res-, 
pirator 71 is connected to the tracheal cannula. A gas mixture of O2 and mixed by a gas mixture 72 and 
dampened by a humidifier 73 is supplied to this artificial respirator 71 . In artificial respiration, the oxygen con- 
centration (RO2) in the inhalation gas to the rat 64 is adjusted, so that the amount of oxygen supplied to the 
rat 64 is controlled. 

In this apparatus arrangement, pulse light incident from the Ti:Sapphire laser 61 on the rat 64 and pulse 
light detected by the streak camera 66 have profiles as shown in the graph of Fig. 13. Time [ns] is plotted along 
the abscissa of this graph, and the detection light Intensity [a.u.] is plotted along the ordinate thereof. A profile 
X is a profile of pulse light incident into the oral cavity of the rat 64. and a profile Y is a profile of detection 
light spreading in the head. An mean flight time <t> between the timing of incidence of the pulse light on the 
rat 64 and the detection timing corresponds to a time <t> between zero time and the barycentral position of 
the profile Y. An average optical pathlength C<t> obtained by multiplying an mean flight time <t> with a speed 
C of light is obtained on the basis of equation (6) described above. This average optical pathlength measure- 
ment was perfomied while the wavelength of the pulse light incident on the rat 64 was changed within the range 
of 740 to 820 nm. The oxygen concentration Fi02 in the Inhalation gas supplied to the rat 64 was changed in 
orders of 100%. 20%, and 15%. 

Fig. 14 Is a graph showing the measurement results. The wavelength [nm] of the incident pulse light is 
plotted along the abscissa, and the average optical pathlength [cm] is plotted along the ordinate. This average 
optical pathlength is calculated on the basis of equation (6) using the value of the speed C of light is 0.022 
cnr^ps. In the graph, a mark O represents a measurement result when the Fi02 concentration Is 100%; A, 20%; 
and □, 1 5%. As shown in this graph, in the short wavelength range, the lower the FiOs concentration becomes, 
and the smaller the amount of oxygen supplied to the rat 64 becomes, the shorter the average optical path- 
length becomes. A short average optical pathlength indicates a high light absorptton in the head portion of the 
rat 64. This coincides with the Hb and Hb02 light absorption profiles shown in the graph of Fig. 1. That is. in 
the graph of Fig. 14. the incident pulse light is greatly absorbed by the deoxygenated hen>oglobin Hb in the 
short wavelength range, thus exhibiting the Hb spectra as shown in Fig. 1. 

A substitution of the measurement result of this average optical pathlength into equation (18) yields the 
ratio of Hb02 concentration to the Hb concentratton. The resultant concentration ratio is substituted Into equa- 
tion (19) to calculate the degree (SO2) of saturation of oxygen In the blood In the head portion of the rat Table 
2 shows the SO2 value obtained from each measurement result represented by each black plot In the graph 
shown in Fig. 14. These plots represent the measurement results with respect to the incident pulse light coay- 
ponents having the wavelengths of 740, 755, and 800 nm. 
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[Table 2] 



Plot 


FiOar/o] 


Average Optical PathLength C<t> [nm] 


SO2 Value 


740 


755 


800 


0 


100 


9.1 


8.7 


8.8 


75 


A 


20 


8.6 


8.1 


8.7 


54 


□ 


15 


8.5 


8.0 


9.1 


39 



As can be understood from the above table, the calculated SO2 value reflects the inhalation oxygen con- 
centration (FIO2). When the concentration Fi02 decreases, the SO2 value also decreases. 

The same effect as in each previous embodiment can be obtained in the third embodiment. The degree 
of saturation of oxygen In the blood of the living tissue can be easily and accurately measured with a simple 
apparatus arrangement 



Claims 

1 . A method of measuring a concentration of an absorptive constituent in a scattering medium, characterized 
by comprising causing light having a predetermined wavelength to be incident from a predetermined pos- 
ition of a scattering medium containing n (n s 2) types of absorptive constituents, detecting the light, scat- 
tering In said scattering medium, at a predetermined position spaced apart from the light incident position 
by a predetermined distance, calculating an average optical pathlength or a mean flight time of the light 
having the predetermined wavelength in said scattering medium on the basis of the detected light, ob- 
taining respective average optical pathlengths or respective mean flight times of light components having 
n+1 types of wavelengths while a wavelength of light Inddent on said scattering medium Is changed, ob- 
taining respective average optical pathlength differences or respective mean flight time differences be- 
tween respective pairs of light components having different wavelengths from the average optical path 
lengths or the mean flight times, and obtaining a concentration ratio of the respective absorptive constit- 
uents in said scattering medium in accordance with a relationship representing that the average optical 
pathlength difference or the mean flight time difference is in Inverse proportion to a difference in absor- 
bance of the absorptive component between the light components having two wavelengths. 

2. A method of measuring a concentration of an absorptive constituent in a scattering medium according to 
claim 1, characterized in that said scattering medium contains two types of absorptive constituents A and 
B, and the incident light components have three different wavelengths AX, k2, and X3, and 

a ratio Va : Vb of a concentration Va and to a concentration Vb of said absorptive constituents A 
and B In said scattering medium is set as follows: 

Va:Vb = {(ebu - eB3ii)A<t>x2,xi " (eew - eBxi)A<t>x3Ai}- { - (caxs - eAxi)A<t>x2.xt + (€ax2 • 

eAXl)A<t>x3;Ll} 

where A<i>J^^ is a mean flight time difference between the incident light components having the 
wavelengths X2 and XI, A<t>x3.M is a mean flight time difference between the incident light components 
having the wavelengths X3 and XI , eaxi and cbxi are extinction coefficients of said absorptive constituents 
Aand B with respect to the incident light component having the wavelength X1 , bau and e^n are extinction 
coefficients of said absorptive constituents A and 8 with respect to the incident light component having 
the wavelength X2, and 8ax3 and sbxs are extinction coefficents of said absorpth^e constituents A and B 
with respect to the incident light component having the wavelength X3. 

3. A method of measuring a concentration of an absorptive constituent in a scattering medium according to 
claim 1 , characterized in that the average optical pathlength is calculated by detecting a light intensity of 
each pulse light obtained at the detection position a predetermined delay time after each pulse light is 
incident on said scattering medium, dividing a sum of products of the delay times and the light intensities 
respectively corresponding to the delay times by a sum of the respective Intensities to obtain a quotient, 
and multiplying the quotient with a speed of light. 

4. A method of measuring a concentration of an absorptive constituent In a scattering medium according to 
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claim 2, characterized in that the average optical pathlength is calculated by detecting a light intensity of 
each pulse light obtained at the detection position a predetermined delay time after each pulse light is 
Incident on said scattering medium, dividing a sum of products of the delay times and the light intensities 
respectively corresponding to the delay times by a sum of the respective intensities to obtain a quotient, 
and multiplying the quotient with a speed of light 

A method of measuring a concentration of an absorptive constituent in a scattering medium according to 
claim 1, characterized in that the average optical pathlength is calculated by causing pulse light to be in- 
cident on said scattering medium, detecting each of photons reaching the detection position with different 
delay times upon incidence of each pulse light, dividing a sum of the delay times of the photons detected 
by a number of photons detected to obtain a quotient, and multiplying the quotient with a speed of light. 

A method of measuring a concentration of an absorptive constituent in a scattering medium according to 
claim 2. characterized in that the average optical pathlength is calculated by causing pulse light to be in- 
cident on said scattering medium, detecting each of photons reaching the detection position with different 
delay times upon incidence of each pulse light, dividing a sum of the delay times of the photons detected 
by a number of photons detected to obtain a quotient, and multiplying the quotient with a speed of light 

A method of measuring a concentration of an absorptive constituent in a scattering medium according to 
claim 1. characterized in that said scattering medium Is a living tissue, said absorptive constituents are 
deoxygenated hemoglobin Hb and oxyhemoglobin HbOa in blood in the living tissue, the incident light has 
a wavelength range of 600 nm to 1.5 ^m, and a degree of saturation of oxygen in the blood in the living 
tissue is obtained on the basis of a calculated ratio of the absorptive constituent concentrations. 

An apparatus for measuring a concentration of an absorptive constituent In a scattering medium, char- 
acterized by comprising a light source for generating pulse light having a predetermined wavelength every 
predetermined time interval, light guide means for causing each pulse light from said light source to be 
incident on a scattering medium containing n (n ^ 2) types of absorptive constituents from a predeter- 
mined position and extracting each pulse light scattering within said scattering medium from a predeter- 
mined position spaced apart from the light incident position by a predetermined distance, photodetecting 
means for time-resolved measuring each pulse light extracted by said light guide means, wavelength con- 
trol means for controlling a wavelength of the pulse light emitted from said light source, optical pathlength 
or flight time calculating means for calculating an average optical pathlength or a mean flight time of light 
components having n+1 wavelengths in said scattering medium in accordance with time-resolved meas- 
urement results of pulse light components having n+1 wavelengths obtained by causing said wavelength 
control means to control the wavelength of the pulse light emitted from said light source, and concentra- 
tion calculating means for obtaining an average optical pathlength difference or a mean flight time dif- 
ference between light components having two wavelengths from the average pathlengths or the mean 
flight times calculated by said optical pathlength or flight time calculating means and obtaining a ratio of 
the absorptive constituent concentrations in said scattering medium in accordance with a relationship rep- 
resenting that the average optical pathlength difference or the mean flight time difference is in Inverse 
proportion to a difference in absorbance of the absorptive component between the light components hav- 
ing two wavelengths. 

An apparatus for measuring a concentration of an absorptive constituent In a scattering medium according 
to claim 8, characterized In that said photodetecting means detects a light Intensity of each pulse light 
obtained at the detection position each predetermined delay time after the pulse light is incident on said 
scattering medium, and said optical pathlength or flight time calculating means divides a sum of products 
of the respective delay times and the light intensities respectively conresponding to the delay times by a 
sum of the light intensities to obtain a quotient, and multiplying the quotient with a speed of light to obtain 
the average optical pathlength. 

An apparatusfor measuring aconcentration of an absorptive constituent in a scattering medium according 
to daim 8, characterized in that said photodetecting means calculates the average optical pathlength by 
detecting each of photons reaching the detection position with different delay times upon incidence of 
each pulse light, and said optical pathlength or flight time calculating means divides a sum of the delay 
times of the photons detected by a number of photons detected to obtain a quotient, and multiplying the 
quotient with a speed of light 
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11. An apparatus for measuring a concentration of an absorptive constituent in a scattering medium according 
to claim 8, characterized in that said light source emits pulse light components having three different wa- 
velengths AX, X2. and X3, and said scattering medium contains two types of absorptive constituents A 
and B, and 

5 said concentration calculating means sets a ratio Va : Vb of a concentration Va to a concentration 

Vb of said absorptive constituents A and B in said scattering medium as follows: 

Va^Vb = {(ebu - eBxi)A<t>x2,xi - (£8X2 - eBxi)A<t>x3,xi} : { - (caxs - eAXl)A<t>x2A1 + (eAX2 - 

eAxi)A<t>x3.xi} 

where A<t>x2;Li is a mean flight time difference between the incident light components having the 
10 wavelengths X2 and X1, A<t>x3,xi is a mean flight time difference between the Incident light components 

having the wavelengths X3 and X1, and baxi a,d ebxi are extinction coefficients of said absorptive constit- 
uents A and B with respect to the incident light component having the wavelength XI, and eAX2 and ^bxz 
are extinction coefficients of said absorptive constituents A and B with respect to the incident light conrv 
ponent having the wavelength X2, and 8ax3 and esxa are extinction coefficients of said absorptive constit- 
15 uents A and B with respect to the incident light component having the wavelength X3. 

1 2. An apparatus for measuring a concentration of an absorptive constituent in a scattering medium according 
to claim 11, characterized in that said light source emits pulse light having a wavelength range of 600 nm 
to 1.5 ^m, said scattering medium is a living tissue, said absorptive constituents are deoxygenated he- 

20 moglobin Hb and oxyhemoglobin Hb02 in blood in the living tissue, and said apparatus further comprises 

SO2 value calculating means for obtaining a degree of saturation of oxygen in the blood of the living tissue 
on the basis of a ratio of said absorptive constituent concentrations output from said concentration cal- 
culating means. 

25 13. Amethod of determining the concentration of an absorptive constituent in a scattering medium, the meth- 
od comprising irradiating the medium at a predetermined position of incidence with radiation having a plur- 
ality of wavelength values, detecting the radiation after passage through the medium at a predetermined 
detection position in respect of each wavelength value and detemnining therefrom a relationship between 
the average optical pathlengths corresponding to the respective wavelengths, thereby to determine the 

9^ concentration of said constituent. 
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Fig . 9 
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